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Abstract 


In  this  report  we  describe  our  analysis  of  the  influence  of  sea  clutter  non 
stationarity  on  the  clutter  covariance  matrix  estimation  and  its  impact  on  the  CFAR 
property  of  the  normalized  adaptive  matched  filter  (NAME).  Three  estimators  have  been 
considered  in  the  analysis,  i.e.  the  sample  covariance  matrix  (SCM),  the  normalized 
sample  covariance  matrix  (NSCM),  and  the  fixed  point  (PF)  estimators.  The  impact  of  non- 
stationarity,  that  emerges  in  the  statistical  analysis  of  the  HH  and  VV  polarized  data,  is 
measured  in  terms  of  differences  between  NAME  nominal  probability  of  false  alarm  (Pfa) 
and  probability  of  detection  (Pd),  and  estimated  ones  measured  by  processing  real  clutter 
data  recorded  by  the  IPIX  radar. 
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1.  Introduction 


To  mitigate  the  deleterious  effeets  of  elutter  and  jammer,  modern  radars  have  adopted 
adaptive  proeessing  teehniques  sueh  as  eonstant  false  alarm  rate  (CFAR)  deteetors, 
adaptive  arrays,  and  spaee  time  adaptive  proeessing  (STAP)  (see  for  instanee  [Kel86], 
[Him98],  [Kra05],  [Mel96],  [Mel97],  [Ran04],  [Ric96],  [War94],  and  referenee  therein).  In 
a  typieal  adaptive  radar  system,  the  disturbance  covariance  matrix  is  estimated  using 
secondary  data,  namely  returns  from  range  cells  spatially  close  to  the  cell  under  test  (CUT) 
and  sharing  the  same  spectral  properties  [Kel86].  In  practice,  these  techniques  are  very 
restrictive  because  they  require  the  environment  to  “remain  stationary  and  homogenous” 
during  adaptation,  which  is  not  always  the  case.  In  fact,  the  secondary  data  are  often 
contaminated  by  interfering  targets,  large  discrete  and  spiky  clutter  and  other  outliers  of 
different  types  making  them  non-homogenous.  The  deleterious  impact  of  this  kind  of  non¬ 
homogeneity  on  the  detection  performances  is  widely  reported  in  the  literature  and  several 
solutions  have  been  proposed  [Ger02],  [Ran05].  However,  another  kind  of  non¬ 
homogeneity  is  caused  by  the  non-stationary  nature  of  the  clutter  which  also  has  a  harmful 
effect  on  the  radar  performances.  Recent  works  on  real  data  revealed  that  the  sea  clutter  is 
a  non  stationary  process  [Hay02],  [Gre04].  Here,  we  propose  to  analyze  this  effect  and  its 
impact  on  CFAR  behavior  of  one  of  the  mostly  used  adaptive  radar  detector,  the 
Normalized  Adaptive  Matched  Filter  (NAMF). 
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2.  Background 

The  problem  of  radar  deteetion  in  non  Gaussian  sea  elutter  has  reeeived  eonsiderable 
attention  in  reeent  past  and  eurrently,  a  substantial  bulk  of  work  is  available  in  the  open 
literature  about  this  topie  [Con94],  [Con95],  [Gin97].  Experimental  data  as  well  as 
physieal  and  theoretieal  arguments  indieate  that  the  elutter  ean  generally  be  modeled  as  a 
eompound-Gaussian  proeess  and  also  show  that  a  satisfaetory  fit  of  the  elutter  amplitude 
probability  density  funetion  (apdf)  at  low  grazing  angles  ean  be  aehieved  through  families 
of  distributions  eontaining  a  shape  parameter  in  addition  to  the  seale  one  [Far97], 
[Con02b],  [Gre07].  Among  these,  the  most  eommonly  adopted  are  the  Weibull  and  K 
distributions,  that  are  partieular  eases  of  the  eompound-Gaussian  family. 

Conventional  radar  deteetors,  namely  those  designed  to  deteet  targets  embedded  in 
Gaussian  elutter,  show  eonsiderable  performanee  degradation  in  the  presenee  of  impulsive 
noise,  even  in  the  ease  of  perfeet  knowledge  of  the  elutter  distribution  parameters  [Con94]. 
Optimized  deteetion  struetures  have  been  proposed  and  assessed  with  referenee  to  eoherent 
pulse  trains  embedded  in  Weibull  or  K  distributed  elutter  with  known  speetral  properties. 
However,  implementation  of  these  deteetion  struetures  requires  knowledge  of  elutter 
statisties  up  to  the  relevant  distribution  parameters,  whieh  is  elearly  unrealistie  in  praetieal 
situations.  Thus,  one  is  interested  in  eanonieal  reeeivers,  namely,  deeision  statisties 
funetionally  independent  of  the  elutter  distribution  parameters  (under  the  noise-only 
hypothesis)  and  whose  deteetion  threshold  is  itself  independent  of  the  elutter  statisties.  A 
first  step  towards  this  task  is  the  reeeiver  introdueed  in  [Con95],  ealled  the  normalized 
matehed  filter  (NMF)  whieh  ensures  the  CFAR  property  with  respeet  to  elutter  distribution 
(in  the  family  of  eompound-Gaussian  models).  This  deteetor  has  also  been  derived 
assuming  Gaussian  disturbanee;  Kraut  et  alii  in  [Kra05]  demonstrated  that  it  is  the  GFRT 
and  Uniformly-Most-Powerful-Invariant  test  for  deteeting  a  target,  know  up  to  a 
multiplieative  faetor,  in  Gaussian  noise  whose  eovarianee  matrix  is  known  but  whose 
power  level  is  unknown.  Henee,  in  order  to  eome  up  with  eompletely  adaptive  deteetion 
struetures  against  a  baekground  of  eompound-Gaussian  elutter,  it  is  neeessary  to  replaee  in 
the  NMF  the  eovarianee  matrix  with  a  suitable  estimate.  The  deteetor  obtained  is  ealled  the 
Normalized  Adaptive  Matehed  Filter  (NAME)  [Ran05]. 

Several  estimators  have  been  proposed  to  estimate  the  eovarianee  matrix  from 
seeondary  data  eolleeted  from  eells  surrounding  the  eell  under  test.  Hereafter,  we  will 
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briefly  review  the  mostly  used  in  practice.  All  these  estimators  are  based  on  the  hypothesis 
that  the  secondary  vectors  do  not  contain  interference  or  targets  and  share  the  same 
covariance  matrix  as  the  primary  data,  the  cell  under  test  (CUT)  under  the  null  hypothesis. 

Before  introducing  the  most  popular  covariance  matrix  estimators,  it  is  useful  to  refer  to 
the  problem  in  its  entirety.  Thus,  the  problem  of  detecting  a  target  signal  in  additive  clutter 
can  be  posed  in  terms  of  the  following  binary  test: 

r//„:z  =  c  z,=c,  i  =  \,...,K 

[//j:z  =  s  +  c  z,=c,  i  =  \,...,K 

where  z,  s  and  c  are  the  A-dimensional  complex  vectors  of  the  samples  from  the  baseband 
equivalents  of  the  received  signal,  the  target  and  the  clutter  in  the  CUT  respectively.  The 
set  of  Zi,  Z2,...,  Zjf  denote  the  A-dimensional  secondary  vectors  assumed  free  of  signals  and 
interferences.  The  useful  signal  s  can  be  modeled  as  s  =  «p  where  p  is  the  target  steering 
vector  and  a  is  an  unknown  parameter  accounting  for  the  channel  propagation  effect  and 
the  target  radar  cross  section.  Under  the  hypothesis  of  compound-Gaussian  clutter  each 
vector  c.  and  c  can  be  modeled  as 


c  =  Vrx , 


where  {x,Xj,...,Xj^}  is  a  sequence  of  independent,  identically  distributed  (IID),  complex, 
zero-mean,  circular  symmetric  Gaussian  vectors  (speckle)  with  unit  power  and  finite, 
positive  definite  covariance  matrix,  in  short  notation  x,x.  e  CA(0,M^)  .  r.  and  r  (texture) 

are  real,  non-negative,  random  variables  independent  of  Xi  and  x. 

Finally,  for  future  convenience,  we  report  here  the  general  structure  of  the  NAMF  test: 


(2) 


where  M  is  any  estimate  of  M  based  upon  the  secondary  vectors,  namely: 
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M=M(Zj,...,z^) 


(3) 


and  M  =  is  the  covariance  matrix  of  the  clutter  normalized  with  respect 

cr  ^  ’  (j 


to  its  power. 
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3.  Covariance  matrix  estimators 

In  this  paragraph  we  describe  three  of  the  most  popular  estimators  of  the  clutter 
covariance  matrix. 

a)  Sample  covariance  matrix:  SCM. 

It  is  an  estimator  of  and  it  is  given  by: 


Mscm  =  ^  Z  ^  )  • 

A  A 


(4) 


Plugging  the  i’^to  (2)  we  obtain  a  detector  that  is  CFAR  only  with  respect  to 

and  whose  probability  of  false  alarm  (Pfa)  strongly  depends  on  the  distribution  of  the 
texture  [Gin02],  [Con02a],  [Con02c].  The  estimator  (4)  is  ML  when  the  clutter  is  Gaussian 
distributed,  then  when  the  texture  values  r.  are  all  equal  and  deterministic.  It  is  also  the 

ML  estimator  when  the  are  random  but  completely  correlated,  i.e.  r  =  rj=...  =  r^ 
[Ric96].  The  performance  of  this  matrix  estimator  in  presence  of  compound-Gaussian 
clutter  has  been  investigated  in  [Gin99]. 

b)  Normalized  sample  covariance  matrix:  NSCM. 

It  is  an  estimator  of  M  and  it  is  given  by: 


M 


NSCM 


A 

K 


K 


z,z,. 


.H 


i=l  * 


A 

K 


i=l 


xfx,. 


i 


(5) 


It  is  very  similar  to  the  SCM  estimate  but  it  uses  normalized  secondary  data,  then  its 
distribution  does  not  depend  on  the  texture  pdf  but  only  on  the  speckle  distribution. 
Plugging  the  (2)  we  obtain  a  detector  that  is  CFAR  with  respect  to  the 

statistics  of  the  texture  and  to  the  clutter  power.  Unfortunately  the  normalization  with 
respect  to  zfz./A  makes  the  NAMF  detector  no  longer  CFAR  with  respect  to  M  [Gin02], 
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[Con02a],  [Con02c].  The  performance  of  the  NSCM  matrix  estimator  in  presence  of 
compound-Gaussian  clutter  has  been  investigated  in  [Gin99]  and  [BauOV]. 

c)  Fixed  point  and  maximum  likelihood  estimator 

The  fixed  point  estimator  (FPE)  arises  as  an  approximation  of  the  ML  estimation 
[Gin02]  (it  has  been  termed  FPE  in  [PasOV]).  It  is  defined  as  a  fixed  point  of  the  function; 


fxM  • 


(D 

A^ 

A 

K 

tr  z"A-' 


z,- 


(6) 


where  (D  =  |  A  e  (C)  |  =  A,  A  positive  defmite|  with  (C)  =  {NxN  matrices  with 

elements  in  C}.  As  shown  in  [PasOV],  equation  M  has  a  solution  of  the  form 

aM ,  where  a  is  an  arbitrary  scaling  factor.  The  only  solution  M  satisfying  the 
normalization  7>^M  'Mj  =  A  is  called  the  Eixed  Point  (EP)  estimate,  . 

In  other  words,  M^p  is  the  unique  solution  of 


=  — 


NF 


L- 


z.zf 


N_j^_ 


x,xf 


AtrzfMpX.  AtTxfMp^x, 


(V) 


satisfying  7>|m  'Mppj  =  A.  Notice  that  Mpp  does  not  depend  on  the  texture  values  T^. 

This  property  is  of  great  importance  in  practice.  As  verified  in  [PasOV],  the  EP  is  the  ML 
estimate  when  the  texture  is  assumed  to  be  an  unknown  deterministic  parameter  and  an 
Approximate  Maximum  Likelihood  (AML)  estimate  when  the  texture  is  assumed  to  be  a 
non  negative  random  variable  [Gin02]. 

To  find  the  EP  solution,  we  must  resolve  the  transcendental  equation  (V).  In  [Gin02]  the 
authors  developed  an  iterative  algorithm  to  calculate  this  solution  in  the  case  of  random 
texture.  The  iterative  algorithm  is  given  by  the  following  equation 


V 


(8) 


-  N  ^ 

Z=1 


z.zf 


where  k  is  the  iteration.  A  deep  analysis  of  this  algorithm  and  its  eonvergenee  rate  as  a 
funetion  of  the  number  of  seeondary  veetors  K,  the  number  of  pulses  N,  and  the  shape  of 
the  elutter  eovarianee  funetion  ean  be  found  in  [Gin02]  and  [PasOV].  Generally  the  number 
of  iterations  Nu  neeessary  to  guarantee  the  eonvergenee  is  small  (say  2  or  3).  In  [Gin02]  it 
was  also  verified  by  simulation  that  the  NAMF  that  uses  the  AML  matrix  estimator  is 
almost  insensitive  to  the  shape  of  the  eovarianee  matrix,  then  it  ean  be  eonsidered  CFAR 
with  respeet  to  the  texture  distribution  and  very  robust  with  respeet  to  the  eovarianee 
matrix.  In  [Con02c]  the  authors  proved  that  if  the  elutter  speetrum  is  symmetrie  about  the 
frequeneyy=0  and  the  initialization  matrix  of  the  iterative  proeedure  (8)  is  done  by 


2^(0)  = 

1^  ^ 


1  K  ,(i),(i)r 


K 


(9) 


where  zf^  and  zf  ^  are  the  real  and  imaginary  part  of  the  eomplex  veetor  z.  respeetively, 
and  is  the  (/,/)th  element  of  the  NxN  matrix  A,  then  the  NAMF  is  CFAR  also  with 

respeet  to  M.  This  is  not  tme  if  the  speetmm  is  not  symmetrieally  shaped  around /=0,  as  it 
is  always  in  sea  elutter  seenarios. 


8 


4.  Statistical  data  analysis 

The  goal  of  this  technical  report  is  to  investigate  the  CFAR  property  of  the  NAMF 
detector  of  eq.  (2)  with  respect  to  the  covariance  structure  and  the  texture  statistics  using 
each  of  the  three  matrix  estimators  reviewed  in  the  previous  paragraph  in  real  sea  clutter 
data,  recorded  by  the  IPIX  radar.  Table  1  shows  the  specifications  of  the  analyzed  files 
19980223_165836_antstep,  19980227_213016_antstep  and  19980227_213808_antstep. 

The  number  of  range  cells  for  the  file  19980223_165836_antstep  is  34,  but 
unfortunately,  the  last  5  (30-34)  contain  a  target  (maybe  a  small  object  floating  on  the 
waves),  then  only  29  cells  can  be  considered  in  our  analysis.  In  the  other  files  there  is  not 
target  and  the  range  cells  are  28. 


Name  of  the 
data  set 

19980223_165836_antstep 

1 9980227_2 1301 6_antstep 

1 9980227_2 1 3808_antstep 

Date,  time  of 
acquisition 

23/02/1998  16:58:36 

27/02/1998  21:30:16 

27/02/1998  21:38:08 

#Range  cells 

34 

28 

28 

Start  range 

3201 m?? 

Range 

resolution 

30  m 

30  m 

15  m 

Pulse  width 

200  ns 

200  ns 

100  ns 

Total  # 
sweep 

60000 

60000 

60000 

Sample  for 
cell 

60000 

Sampled  at  30  m 

60000 

Sampled  at  30  m 

60000 

Sampled  at  15  m 

PRF 

1  KHz 

1  KHz 

1  KHz 

Frequency 

RF 

9.39  GHz 

9.39  GHz 

9.39  GHz 

Radar  and 

wave 

geometry 

N 

\ 

N 

N 

S 

■ 

S 

Table  1  -  Characteristics  of  the  analyzed  files. 
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The  data  have  been  previously  proeessed  in  order  to  remove  the  dc  offset  and  the  phase 
imbalance  due  to  hardware  imperfections  and  then  stored  in  an  N,  x  complex  matrix.  A 

detailed  statistical  analysis  of  adopted  real  data  has  been  performed  as  in  [Gre06]  and  in 
[Gin06].  The  results  have  shown  a  reasonable  fit  of  the  data  to  the  K  model. 

The  shape  parameter  v  is  not  constant  on  all  the  range  cells.  In  Figs.  1-3  we  show  the 
estimated  values  for  VV  and  HH  data  for  the  three  files.  As  expected  the  values  of  v  for 
the  HH  data  are  always  lower  than  that  for  the  W  data,  then  the  HH  polarization  is  spikier 
than  the  VV.  The  ratio  between  the  v  parameters  of  the  W  data  and  the  HH  data  is  close 
to  0.9  for  each  cell  and  each  analyzed  file.  The  mean  value  of  v  for  the  VV  data  and  HH 
data  are  reported  in  table  2.  The  spikiest  data  are  those  of  file  19980227_213016_antstep. 


File  name 

W 

HH 

19980223_165836_antstep 

0.7563 

0.4355 

19980227_213016_antstep 

0.6434 

0.3981 

1 9980227_2 1 3  808_antstep 

1.0046 

0.5141 

Table  2  -  Vniean  for  the  analyzed  files. 


Fig.  1  -  File  19980223_165836_antstep. 
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shape  parameter  v  shape  parameter  v 


Fig.  2  -  File  19980227_213016_antstep. 


Fig.  3  -  File  19980227_213808_antstep. 
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5.  Spectral  analysis 

We  estimated  also  the  powers  speetral  density  (PSD)  of  the  clutter  in  each  range  cell  by 
means  of  Welch  modified  periodogram  with  rectangular  windows  of  256  samples  and  an 
overlap  of  50  %.  The  results  are  reported  in  Figs.  4  and  5  for  the  D'  file 
(19980223_165836_antstep),  where  the  periodogram  normalized  with  respect  to  the 
estimated  power  for  the  29  cells  is  plotted  in  a  color  scale  from  blue  to  red,  as  a  function  of 
the  number  of  cell  and  the  normalized  (with  respect  to  the  radar  PRF)  frequency.  The 
spatial  non  stationarity  of  the  clutter  is  evident,  the  PSD  varies  from  cell  to  cell.  The  peak 
of  the  PSD  is  generally  between /=0. 05  and/=0.1,  as  verified  also  in  Fig.  6  where  the 
average  PSD  is  reported  for  both  polarizations.  The  average  PSD  has  been  calculated  as 
the  mean  of  the  periodograms  estimated  on  the  29  cells. 

We  obtain  similar  results  for  the  2"‘^  and  the  3’^‘^  files  as  shown  in  Figs.  7-9  for  file 
19980227_213016_antstep  and  Figs. 10-12  for  file  19980227_213808_antstep. 


range  cell  number 

Fig.  4  -  PSD  of  sea  clutter  as  a  function  of  range  cell  number,  VV  data,  P‘  file. 
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normalized  frequency 


5  10  15  20  25 


range  cell  number 

Fig.  5  -  PSD  of  sea  clutter  as  a  function  of  range  cell  number,  HH  data,  file. 


normalized  frequency 


Fig.  6  -  Average  PSD  of  sea  clutter,  file. 
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normalized  frequency 


5  10  15  20  25 


range  cell  number 

Fig.  7  -  PSD  of  sea  clutter  as  a  function  of  range  cell  number,  VV  data,  2'“*  file. 
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Fig.  8  -  PSD  of  sea  clutter  as  a  function  of  range  cell  number,  HH  data,  2'“*  file. 
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normalized  frequency 


normalized  frequency 


Fig.  9  -  Average  PSD  of  sea  elutter,  2"‘*  file. 
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Fig.  10  -  PSD  of  sea  elutter  as  a  funetion  of  range  eell  number,  W  data,  3^'*  file. 
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normalized  frequency 
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Fig.  11  -  PSD  of  sea  clutter  as  a  function  of  range  cell  number,  HH  data,  3'^'*  file. 


normalized  frequency 


Fig.  12  -  Average  PSD  of  sea  clutter,  3^'*  file. 
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To  complete  the  spectral  analysis,  we  ealeulated  also  the  speetrogram  of  the  data  for 
eaeh  range  cell  as  in  [Gin06].  For  the  l'^*  file  we  show  here  only  the  results  for  the  1 1*’’  eell, 
in  Figs.  13  and  14  for  VV  and  HH  data  respectively.  In  both  the  polarizations  the  behavior 
is  similar  and  some  temporal  periodicities  are  apparent,  particularly  in  the  ranges  5-20  sec 
and  35-50  sec.  Again  the  temporal  non-stationary  behavior  of  the  sea  clutter  is  evident. 

For  the  file  we  report  in  this  section  the  results  for  the  4*  and  11*  range  cell  in 
Figs.  15-20.  For  the  3'^‘^  file  we  report  the  7*,  13*  and  27*  in  Figs.  21-25  for  both 
polarizations.  Also  in  these  two  files  the  temporal  non-stationarity  is  apparent. 


time(sec) 


Fig.  13  -  Spectrogram  of  sea  clutter,  1 1*  cell,  VV  data,  file. 
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normalized  frequency  normalized  frequency 


time(sec) 

Fig.  14  -  Spectrogram  of  sea  clutter,  1 1*  cell,  HH  data,  1®“  file. 

#1 


Fig.  15  -  Spectrogram  of  the  1®‘  cell,  W  data,  2"^  file. 
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#1 


0.22 
0.2 

0.18 

0.16 

0.14 

0.12 

0.1 

0.08 

0.06 

0.04 
0.02 

0  10  20  30  40  50 

time(sec) 


Fig.  17  -  Spectrogram  of  the  4*  cell,  W  data,  2"“^  file. 
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normalized  frequency  normalized  frequency 
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Fig.  18  -  Spectrogram  of  the  4*  cell,  HH  data,  2"“^  file. 
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Fig.  19  -  Spectrogram  of  the  1 1*  cell,  W  data,  2"“^  file. 
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normalized  frequency 
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Fig.  20  -  Spectrogram  of  the  1 cell,  HH  data,  2"“^  file. 
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Fig.  21  -  Spectrogram  of  the  7*  cell,  W  data,  file. 
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normalized  frequency 
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Fig.  22  -  Spectrogram  of  the  7*  cell,  HH  data,  file. 
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Fig.  23  -  Spectrogram  of  the  13*  cell,  W  data,  3'“*  file. 
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normalized  frequency 
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Fig.  24  -  Spectrogram  of  the  13*  cell,  HH  data,  S"'*  file. 
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Fig.  25  -  Spectrogram  of  the  27*  cell,  W  data,  S"'*  file. 
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#27 


time(sec) 

Fig.  26  -  Spectrogram  of  the  27*  cell,  HH  data,  file. 


Summarizing  our  results  on  the  statistical  analysis  on  sea  clutter  file  we  can  state  that: 

1 .  The  sea  clutter  show  a  good  fit  to  K  model  but  the  shape  parameter  of  the  distribution 
changes  from  cell  to  cell. 

2.  The  speckle  PSD  is  not  constant  in  time  and  space,  so  the  clutter  is  not  spatially  and 
temporally  stationary.  The  spectrogram  evidences  some  temporal  periodicity  in 
spectrum  PSD  behavior. 

The  aim  of  the  analysis  described  in  the  following  sections  is  to  verify  and  measure  the 
impact  of  these  spatially  and  temporally  non-stationarity  on  the  CFARness  and  in  general 
on  the  performance  of  the  NAMF. 
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6.  Performance  analysis  and  CFAR  property 

The  procedure  used  to  assess  the  performances  is  now  described.  We  consider  an 
X  (if  + 1)  data  window,  where  N  is  the  number  of  pulses  and  K  the  number  of  secondary 

cells,  the  primary  cell  (CUT)  is  set  in  the  middle  of  the  window.  The  data  window  is  slid  in 
space  from  range  bin  to  range  bin  and  in  time  with  an  overlap  of  50  %  {N/2  samples)  until 
the  end  of  the  dataset.  The  overall  number  of  trials  is: 


N. 


iN,-N/2) 


trials 


Nl2 


where  W=60000  and  Nc=29.  Moreover,  for  evaluating  the  performance  of  the  NAMF,  we 
set  N=S  and  K=\6.  The  steering  vector  of  the  signal  is: 


P  = 


1,  exp  {j27rf^ ) , . . . ,  exp  (72;z-  ( A  - 1)  /^  )]^ 


where  is  the  normalized  Doppler  frequency  of  the  target. 

In  analyzing  the  behavior  of  the  NAMF  with  the  real  data  our  aim  is  to  highlight  the 
impact  of  clutter  non-stationarity  on  the  false  alarm  rate  (FAR).  As  written  in  Section  3,  we 
know  that: 

(/)  the  NAMF  with  the  SCM  is  CFAR  with  respect  to  Me  and  to  M,  but  not  with  respect 
to  the  statistics  of  the  texture; 

(n)  the  NAMF  with  the  NSCM  is  CFAR  with  respect  to  the  statistics  of  the  texture,  but 
not  with  respect  to  M; 

(in)  the  NAMF  with  the  FP  is  CFAR  with  respect  to  the  statistics  of  the  texture,  and 
very  robust  with  respect  to  M. 

Then,  to  set  the  threshold  for  a  nominal  probability  of  false  alarm  F),^q=10  '  and 
^  in  the  NAMF  with  each  of  the  three  matrix  estimators,  we  generated  K- 
distributed  clutter  with  covariance  matrix  equal  to  the  average  covariance  matrix  of  the  real 
data  and  with  shape  parameter  equal  to  the  estimated  .  The  average  covariance  matrix 
has  been  built  based  on  the  estimated  correlation  as: 
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(10) 


1  N^. 

7777  2  S  7  («)7  («  +  '”) 

^^=1  i=0 

where  z^(«)  in  the  sample  of  the  range  cell  on  the  clutter  complex  envelope/ 

Running  Monte  Carlo  simulation  with  this  simulated  data  for  different  values  of  the  target 
Doppler  frequency  fn,  we  observed  that,  even  if  the  NAMF-NSCM  and  the  NAMF-FP  are 
not  CFAR  with  respect  to  M,  the  nominal  threshold  shows  very  small  variations  as 
function  of  fo  and  it  can  be  practically  considered  constant.  The  nominal  thresholds  are 
reported  in  Table  3  for  the  file. 


PpA-m' 

SCM 

0.537 

0.760 

NSCM 

0.460 

0.682 

FP 

0.485 

0.712 

Table  3a  -  NAMF  threshold  values,  VV  data 


Pfa-1^'' 

Pfa^I^-"- 

SCM 

0.604 

0.816 

NSCM 

0.456 

0.680 

FP 

0.484 

0.710 

Table  3b  -  NAMF  threshold  values,  A=8,  i:=16,  HH  data 


After  setting  the  nominal  thresholds  we  estimated  by  Monte  Carlo  simulation  the 
probability  of  false  alarm  of  the  NAMF  fed  by  the  real  data  in  both  the  polarizations.  The 
results  are  shown  for  the  C*  file  in  Figs.  27-28  for  the  VV  data  and  Figs.  29-30  for  the  HH 
data  as  a  function  of  the  target  Doppler  frequency  fo. 

’  The  mean  value  for  each  range  cell  has  been  subtracted  before  processing  the  data. 
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probability  of  false  alarm  2  probability  of  false  alarm 
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target  frequency  f^ 

27  -  Probability  of  false  alarm  with  real  W  data,  A^=8,  if=16,  P/r^o=10'',  file. 


target  frequency  f^ 


Fig.  28  -  Probability  of  false  alarm  with  real  W  data,  N^8,  K^16,  Pi?^o=10'^,  P‘  file. 
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probability  of  false  alarm  2  probability  of  false  alarm 


target  frequency  f^ 


29  -  Probability  of  false  alarm  with  real  HH  data,  N=8,  K^16,  Pf,4o=10'',  T*  file. 


target  frequency  f^ 


Fig.  30  -  Probability  of  false  alarm  with  real  HH  data,  N=8,  K^16,  P/r^o=10'^,  file. 
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2 

The  differences  between  nominal  and  actual  Pfa  are  evident,  particularly  for  Pfao=^^'  ■ 
The  actual  Pfa  of  NAMF-NSCM  and  NAMF-FP  are  close  to  the  nominal  one  only  for  fo  in 
the  noise  floor  of  the  clutter.  When  fo  is  close  to  the  peak  of  the  PSD,  the  real  Pfa  is  higher 
than  the  nominal  one,  particularly  for  Pfao=^^  ■  On  the  contrary,  the  actual  Pfa  of  the 
NAMF-SCM  is  almost  always  lower  than  the  nominal  one.  The  results  are  similar  for  both 
polarizations.  It  is  important  to  observe  that  the  deviations  from  the  nominal  value  of  the 
Pfa  are  higher  in  the  area  of  the  clutter  spectrum  in  which  the  variations  due  to  the  non- 
stationarity  of  the  clutter  are  greater.  The  variations  of  the  noise  floor  are  negligible.  It  is 
the  peak  of  the  PSD  or  the  Doppler  centroid  that  moves  with  the  long  waves  originating  the 
almost  periodic  behavior  of  the  spectrogram  as  in  Figs.  5-6. 

The  differences  in  the  NAMF-SCM  can  be  mostly  due  to  the  non-stationarity  of  the  shape 
parameter  of  the  K  distribution  to  which  the  data  belong,  more  than  to  the  non-stationarity 
of  the  covariance  matrix.  The  SCM  is  particularly  sensitive  to  the  clutter  texture 
probability  density  function. 

Obtaining  a  Pfa  that  is  lower  than  the  nominal  one  does  not  mean  in  general  having  a 
good  result.  The  performance  of  the  system  depends  also  on  the  probability  of  detection. 
To  evaluate  it,  we  synthetically  generated  a  Swerling  I  target,  so  that  the  data  samples  in 
the  CUT  have  the  following  expression: 

z{n)  =  aQ-K^[j2nfj^n)  +  d{n)  (12) 


for  n  =  0,1,...  A-1 ;  a  is  a  complex  Gaussian  random  variable,  i.e.  «  e  CA'(0,cr^^,yD  is 

the  Doppler  frequency  of  the  target  and  d{n)  is  the  clutter.  The  clutter-to-noise  ratio  (CNR) 
is  defined  as  SCR  =  cr^  /cr^  .  In  the  simulation  with  synthetic  clutter  we  set  crj  =  1 ,  in  the 
test  with  real  data  the  average  power  of  the  clutter  has  been  estimated  from  the  data  as 


N,.  N,,-\ 


^1  = 


NN. 


k=x  ,=0 


(13) 


where  Nc  is  the  number  of  range  cells  and  Nd  the  number  of  data  per  cell. 
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The  results  are  summarized  in  the  Figs.  31-34  for  the  1®^  file,  where  we  set ^=0.  It  is 
evident  from  these  figures  that  the  performanees  of  the  three  deteetors  are  poorer  with  the 
real  non-stationary  data  than  with  the  simulated  stationary  elutter.  The  non-stationarity  of 
the  elutter  infiuenees  also  the  probability  of  deteetion  of  the  NAMF. 

The  performances  of  the  NAMF  for  the  2"‘^  file  are  summarized  in  Figs.  35-38.  For 
a  nominal  ^  the  estimated  Pfa  of  NAMF-NSCM  and  NAMF-FP  are  pretty  close  to 

the  theoretical  one.  This  is  not  true  for  the  NAMF-SCM  that  exhibits  a  Pfa  lower  than  the 
nominal  one.  For  Pfao^^^  all  the  estimated  Pfa  are  different  from  the  nominal  one.  The 
results  are  similar  for  both  polarizations.  From  Figs.  39-42,  it  appear  evident  that  the  actual 
Pd  is  always  lower  that  the  nominal  Pd  for  each  detector  and  matrix  estimator. 


Fig.  31  -  Probability  of  detection,  K^\6,  W  polarization,  file. 
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Fig.  32  -  Probability  of  detection,  N=8,  K^16,  W  polarization,  T*  file. 


Fig.  33  -  Probability  of  detection,  N=S,  K^16,  PpAo^lO'^,  HH  polarization,  1"‘  file. 
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Fig.  34  -  Probability  of  detection,  A^=8,  K^\6,  HH  polarization,  file. 
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Fig.  35  -  Probability  of  false  alarm  with  real  W  data,  N=%,  K=\6,  Piv4o=10'* 


32 


probability  of  false  alarm  probability  of  false  alarm 
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Fig.  36  -  Probability  of  false  alarm  with  real  W  data,  N=8,  K^16, 


target  frequency  f^ 


Fig.  37  -  Probability  of  false  alarm  with  real  HH  data,  N=8,  K^16,  Piv4o=10'* 
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nominal  probability  of  detection  probability  of  false  alarm 
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Fig.  38  -  Probability  of  false  alarm  with  real  HH  data,  N=8,  K^16,  Piv(o=10'^ 


Fig.  39  -  Probability  of  detection,  N=8,  K^16,  P/r^o=10  ',  W  polarization. 
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nominal  probability  of  detection  nominal  probability  of  detection 


Fig.  40  -  Probability  of  detection,  N=8,  K^16,  Pf.4o=10'^,  VV  polarization. 


Fig.  41  -  Probability  of  detection,  N=8,  K^16,  P/r^o=10'*,  HH  polarization. 
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Fig.  42  -  Probability  of  detection,  N=8,  K^16,  Pf.4o=10'^,  HH  polarization. 


We  continued  our  analysis  testing  the  NAMF  19980227_213808_antstep  file  data. 
In  this  file  the  range  resolution  is  15  m.  The  Pfa  and  the  Pd  are  plotted  in  Figs.  43-50. 
Surprisingly,  the  probability  of  detection  for  the  VV  data  is  pretty  close  to  the  nominal  one. 
The  only  difference  with  the  other  data  is  that  this  VV  polarized  clutter  is  less  spiky  than 
the  others  in  the  other  two  files.  We  can  guess  that  the  non-stationarity  have  a  greater 
weight  with  spikier  data.  As  a  matter  of  fact,  the  Pd  on  the  HH  data  shows  same  behaviour 
than  in  the  previously  analyzed  files. 
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probability  of  false  alarm  probability  of  false  alarm 
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Fig.  43  -  Probability  of  false  alarm  with  real  W  data,  N=S,  if=16,  Piv4o=10'* 


target  frequency  f^ 


Fig.  44  -  Probability  of  false  alarm  with  real  W  data,  K^16,  P/7^o=10'^ 
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Fig.  45  -  Probability  of  false  alarm  with  real  HH  data,  N=8,  K^16,  Piv4o=10'* 


target  frequency  f^ 


Fig.  46  -  Probability  of  false  alarm  with  real  HH  data,  N=8,  K^16,  P/7^o=10'^ 
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nominal  probability  of  detection  nominal  probability  of  detection 


Fig.  47  -  Probability  of  detection,  N=8,  K^16,  Pf.4o=10"',  VV  polarization. 


Fig.  48  -  Probability  of  detection,  N=8,  K^16,  Pf.4o=10'^,  W  polarization. 
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nominal  probability  of  detection  nominal  probability  of  detection 


Fig.  49  -  Probability  of  detection,  N=8,  K^16,  Pf.4o=10"',  HH  polarization. 


Fig.  50  -  Probability  of  detection,  N=8,  K^16,  Pf.4o=10'^,  HH  polarization. 
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7.  Conclusions 

In  this  technical  report  we  have  analyzed  HH  and  VY  polarized  lake  clutter  data  and  we 
have  used  them  to  test  the  performance  of  the  Normalized  Adaptive  Matched  Filter 
(NAMF)  detector  with  the  goal  of  verifying  the  impact  of  temporal  and  spatial  non- 
stationarity  on  the  detector  False  Alarm  Rate  (FAR)  and  in  general  on  the  performance  of 
the  NAMF  adopting  different  clutter  covariance  estimators.  We  applied  the  sample 
covariance  matrix  (SCM),  the  normalized  sample  covariance  matrix  (NSCM),  and  the 
fixed  point  (FP)  matrix  estimators.  All  the  known  properties  of  NAMF-SCM,  NAMF- 
NSCM  and  NAMF-FP  are  based  on  the  assumption  that  the  clutter  in  the  primary  and 
secondary  vectors  shares  same  statistical  properties,  that  is,  the  clutter  is  stationary  in  time 
and  space.  In  real  sea  clutter  scenario,  unfortunately,  often  this  is  not  the  case. 
Our  statistical  analysis  shows  that: 

1)  The  sea  clutter  show  a  good  fit  to  K  model  but  the  shape  parameter  of  the 
distribution  changes  from  cell  to  cell. 

2)  The  speckle  PSD  is  no-constant  in  time  and  space,  so  the  clutter  is  not  spatially  and 
temporally  stationary.  The  spectrogram  evidences  some  temporal  periodicity  in 
spectrum  PSD  behavior. 

To  measure  the  impact  of  clutter  non-stationarity  on  CFAR  property  of  the  NAMF,  we 
generated  a  K  clutter  with  a  covariance  matrix  equal  to  the  average  covariance  matrix  of 
the  real  data  and  with  a  shape  parameter  equal  to  the  estimated  .  We  performed  a 

thorough  numerical  analysis  by  processing  these  simulated  data,  for  different  values  of  the 
target  Doppler  frequency  fo,  and  we  set  the  threshold  for  a  nominal  probability  of  false 
alarm  =  10  '  and  =  10^^ . 

After  setting  the  nominal  thresholds  we  estimated  the  probability  of  false  alarm  and  of 
detection  of  the  NAMF  fed  by  the  real  data.  If  the  data  had  same  statistical  properties  of 
the  simulated  data,  the  performance  of  the  NAMF  with  real  data  and  with  simulated  data 
would  be  very  close.  Conversely,  in  our  results  we  observed  large  differences  in  the  Pfa 
and  in  the  Pd,  almost  always.  Most  of  these  differences  can  be  attributed  to  the  spatial  and 
temporal  non-stationarity  of  real  data  that  can  imply  a  variation  of  the  Pfa  of  one  order  of 
magnitude  with  respect  to  the  nominal  one,  and  of  the  SCR  for  a  fixed  Pd  of  2-5  dB.  The 
impact  of  the  nonstationarity  seems  to  be  stronger  and  stronger  with  increasing  spikiness 
and  decreasing  Pfa- 
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